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Abstract 

Kij'envahie  weighting  niipeiirH  in  Homo  noise  s\il).s|>ac,c 
methods,  in  parametric  signal  snhspace  fitting  methods, 
and  in  nonparamctric  snhspnee  adaptive  nnlling  hcamform- 
ing.  For  a  two-source  array  processing  scenario,  normalized 
cigenvaltics’  expressions  Aj  and  Aj  are  reduced  to  forms  de¬ 
pending  only  on  a  r<!al  triplet:  phase-dependent  variable 
phase-independent  variable  17,  and  power  ratio  (^,»/)  is 
confined  to  an  isosceles-like  region.  We  characterize 

•  this  isosceles-like  region  and  the  many-to-one  mapping 
from  the  Cartesian  product  of  the  temporal  and  spatial 
correlation  unit-disks  onto  this  region,  and 

•  the  behavior  of  the  eigenvalues  and  (heir  ratio  as  func¬ 
tions  of  the  real  trijdet  both  analytically  and  graphi- 
rnlly  %vilh  respect  (o  array  [Processing. 

1 .  Tiit.roduclion 


:a  D'iifi'i-finoo 

We  give  sjiecial  treatment  for  eipial  jpow'er  arrivals  be- 
cau.se  of  its  many  unifpie  charac(eristirs  and  its  importance 
in  real  applications. 

'I’hc  analytical  treatment  begins  by  reviewing  the  ex¬ 
pressions  of  the  noise-free  eigenvalues  of  the  quadratic  char¬ 
acteristic  equations  of  the  non-llcrmitian  product  of  the 
temporal  and  spatial  correlation  matrix  given  in  Hudson’s 
text[ld,  [pp.  5Iil-55].  The  eigenvalues  are  normalized  with 
respect  to  the  product  of  the  number  of  sensors  and  the  sen¬ 
sor  level  power  of  the  weaker  source  rro,  i.e.,  tti  >  ttt.  The 
normalized  large  and  small  eigenvalues  expressions  Aj  and 
A2  are  rcducerl  to  forms  depending  only  on  the  real  triplet 
(^,1/,  ^).  Here,  ^  is  the  [Power  ratio  between  the  strong 
to  woat  sources  at  the  sensor  level.  'Flic  real  pair  (^,tj)  arc 
defined  in  terms  of  the  normalized  tenqporal  and  spatial  co¬ 
efficients  p  and  rcs[Pcctively  with  each  constrained  to  a 
unit-disk  through 

v-(i  -  kpr’)(i  -IH-), 


In  this  [pa[per,  we  discuss  some  stre.ss  measure  to  array 
processing  algorithms  a|pphed  toscenaricps  with  two  sources. 
'I'he  word  "stress”  for  a  sci'iiario  is  used  to  indicate  lh(>  diffi¬ 
culty  that  an  algorithm  luis  in  determining  the  direction  of 
arrival  of  each  of  t  he  two  sources.  Schmidt’s  dissertation  [1] 
qualitatively  [pointed  out  thri'e  stressing  factors  to  an  array 
[Processing  algorit  hm: 

1.  For  small  or  medium-sized  sensor  arrays,  the  resolu¬ 
tion  is  not  sufficiently  high  so  we  can  have  unresolved 
arrivals. 

2.  Furthermore,  high  or  100%  source  correlation  such  as 
from  multipaths  will  cause  ill-conditioning  or  rank  de¬ 
ficiency  of  (he  noise-free  data  covariance  matrix. 

3.  The  other  stress  factor  to  a  direction-finding  algorithm 
or  scenario  is  the  strength  rat  io  between  the  two  signals. 

Here,  we  use  eigenvalue  ratios  to  quantitatively  combine  the 
three  strc.ssing  factors.  Eigenvalues  in  weighting  expressions 
appear  in  three  categories  of  array  processing  algorithms: 

•  in  noise  subs[Pacc  methods  [2,3,'1], 

•  in  parametric  signal  subspace  fitting  methods  [<1,5,6, 7, 
8,9,10]  or  its  associated  nonparametric  subspace  beam- 
forming  [1 1],  and 

•  in  nonparametric  subspace  adaptive  nulling  beamform¬ 
ing  [12,13]. 

Only  [2]  is  al/out  weighted  noise-eigenvectors  with  weights 
involving  the  corresponding  noise  eigenvalues.  See  the  re¬ 
cent  result  on  this  in  [3].  All  others  use  weighted  signal 
eigenvectors  with  weights  involving  the  corresponding  sig¬ 
nal  eigenvalues. 


^  =  He(/)Ci')  =  |/■p|[d'|c(ps(argp  -  arge’p). 

Using  (Pur  notation,  the  normalized  eigenvalues  assume  (he 
following  form 


Before  finishing  the  introduction  of  symbols  and  defini¬ 
tions,  we  point  out  the  importance  of  the  sensor  coordinate 
origin.  When  two  arrivals  are  100%.  correlated,  the  phase 
difference  bp'tween  the  two-source  signals  plays  an  impor¬ 
tant  role.  The  specification  of  this  variable  depends  on  the 
choice  of  the  coordinate  origin.  Frequently,  performance 
such  as  the  Cramer-Rao  Lower  Bound  (CRLB)  is  plotted 
against  this  variable.  If  two  different  origins  are  used,  the 
same  results  can  be  seen  differently  but  are  equivalent  with 
a  horizontal  shift  and  possibly  with  the  help  of  27r-modulo 
wra(p-around.  Naturally,  the  covariance  matrix  and  their 
eigenvalues  are  not  dependent  on  the  choice  of  the  origin, 
|put  the  expressions  for  the  phase  angles  of  both  the  steering 
vectors  and  complex  source  amplitudes  of  the  signals  are. 

In  the  special  ca.se  of  plane  waves,  an  array  is  called 
pairwi.se  symmetric  o.  inverse  symmetric  [15]  if  for  each 
sensor  located  away  from  the  origin  at  position  coordinates 
there  is  al.so  a  sensor  located  at  (— r,-,  — y,-,  — z,). 
The  inner  product  of  two  steering  vectors  for  plane  waves 
impinging  onto  pairwise  symmetric  arrays  is  real.  For  such 
case,  the  usual  two-dimensional  unit-disk  associated  with 
spatial  correlation  4>  degenerates  into  a  one-dimensional  one. 
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'lh(iii,  t.lic  <liir<'r<!iu'(:  l)<;I.W('oii  (lio  (.('iiiporal  and  spa¬ 

tial  correlation  cocllicients  is  reduced  to  the  temporal  phase 
dilTereiice  between  the  two  arrivals  measured  at  the  origin 
oT  the  symmetric  array.  While  the  discn.ssions  of  eigenval¬ 
ues  in  this  report  are  applicable  to  generic  sttrering  vectors, 
we  call  1)  the  phase-independent  variable  and  ^  the  phase- 
dependent  variable  for  convenience  even  though  such  names 
were  motivated  by  the  special  case  of  plane  waves  impinging 
onto  pairwise  symmetric  arrays. 

2.  An  Isosceles  Right-Trianglc-Like  Region 

The  phase-independent  variable  r;  depends  only  on  the 
magnitude  of  p  and  <f>.  The  phase-dependent  variable  ^  <le- 
pends  also  on  the  angular  positions  of  p  and  <f>.  The  mapping 
from  the  complex  pair  (p,  4>)  to  the  real  pair  (^,  ij)  is  many  to 
one.  The  range  of  this  mapping  is  an  isosceles  right-triangle¬ 
like  region  bounded  oti  its  left  and  right  liy  two  symmetric, 
parabohis,  r;  =;  (^  ±  1)‘,  and  down  below  by  a  straight  base¬ 
line,  ;;  =  0.  We  note  the  zero  slope  of  each  parabola  where 
it  meets  the  straight  line  and  the  90®  angle  that  the  two 
parabolas  intersect  each  other.  We  characterize  this  region 
with  respect  to  the  temporal  and  spatial  correlation  coefTi- 
cients  of  array  processing  scenarios  in  Figures  1  and  2. 

We  will  see  later  that  among  the  three  parameters 
appearing  in  the  eigenvalues’  expressions,  t; 
is  the  most  important  one.  In  the  following,  we  u.sc  Fig¬ 
ure  1  to  characterize  (he  many-to-one  mapping  from  the 
Cartesian  product  of  the  two  unit-disks  onto  the  isosceles- 
like  region  with  respect  to  rj.  We  note  that  for  given  ij, 
max(|p|,|0|)  =  y^l  —  y?/.  Therefore,  the  apex  point  (0,1) 
corresponding  to  r;  =  1  can  only  come  from  the  two  centers 
of  the  unit-disks,  i.e.,  p  =  =  0.  As  ?;  d<;cr(!)Lses  from  1 

to  0,  i.e.,  we  move  from  the  apex  towards  the  baseline,  we 
shade  the  part  of  the  Cartesian  product  of  the  two  unit- 
disks  mapped  into  the  correspondingly  shaded  area  of  the 
iso.sceles-like  region.  I'herefore,  as  the  shaded  area  of  the 
iso.sceles-like  region  expands  from  the  apex  towards  the  base¬ 
line,  the  shaded  areas  of  the  two  unit-disks  expand  from  the 
centers  toward  the  unit  circumferences. 

When  we  move  one-fourth  of  the  way  toward  the  base¬ 
line,  i.e.,  T]  >  0.75,  we  find  max(|p|,  |<^|)  =  \/l  —  \/b.75  rj 
0.36.  For  two  uncorrelated  plane  wave  arrivals  impinging  on 
a  uniform  linear  array,  the  first  sidelobe  peak  has  a  height 
of  101  «  v/0.05  ss  0.22.  At  T]  =  0.75,  the  two  interferers 
are  already  within  a  beamwidth  of  each  other.  Let  us  shade 
the  corresponding  range  and  domain  of  this  many-to-one 
mapping  using  db®  hatching  lines  for  this  rj  >  0.75  region. 

Similarly,  when  we  have  moved  halfway  toward  the 
baseline  in  the  isosceles-like  region  r]  >  0.50,  we  find 
max(|p|,|0|)  =  \/l  —  w  VO.3  «  0.54.  Let  us  over¬ 
lay  onto  the  earlier  picture  by  shading  the  corresponding 
range  and  domain  using  offset  45°  hatching  lines  for  this 
r]  >  0.50  value.  We  note  that  the  pieviously  45®  hatched 
region  are  double-hatched  now.  For  any  point  in  the  in¬ 
cremental  area  of  the  isosceles-like  region  corresponding  to 
the  singly  45®  hatched  but  not  double-hatched  area,  i.e., 
0.75  >  rj  >  0.50,  at  least  one  of  the  originating  p  or  0  must 
be  located  in  the  correspondingly  shaded  incremental  annu¬ 
lus,  i.e.,  0.36  <  |p|  or  |01  <  0.54. 


Wlu-ii  we  move  three-fourt  h  ol  thi-  way  toward  the  b.i.se- 
liiic,  i.e.,  rj  >  0.25,  we  find  max(|p|,|0|)  =  x/ 1  —  « 

0.7.  The  two  interferers  are  now  within  half  of  the  half- 
power  beamwidt  h  of  each  ot  her. 

We  id.so  iioti'  (  hat  ?/  :=  ( I  —  |0|‘)(  1  —  |p|")  is  the  height  of 
the  unshaded  trapezoidal-like  region.  We  can  associate  rr-rj 
as  the  product  of  (  he  areas  of  the  two  unshaded  annulus  re¬ 
gions,  ir(l  — 10|-)  and  Tr(l  — |p|*).  As  this  height  reduces  from 
I  to  0,  i.e.,  the  unshailed  t  rapezoid  decren.ses  from  the  origi¬ 
nal  whole  isosceles  all  the  way  to  the  zero  thickue.ss  baseline, 
at  least  one  of  the  corresponding  white  annuli  shrinks  from 
the  original  unit-disk  to  a  zero  width  ring. 

When  the  scenario  is  stressful,  rj  is  small,  the  corre- 
spoixling  annuli  with  unit  outer  radius  arc  thin.  The  area 
of  the  temporal  correlation  p  annulus  can  be  approximated 
by 

7r(l  -|p|2)  =  tr(l-P|p|)(l-|p|)«2rr(l-|p|). 

Suppose  the  temporal  correlation  coefTicient  has  a  magni¬ 
tude  of  99%,  then  the  thickness  of  the  p  annulus  is  0.01. 
The  area  of  the  annulus  can  be  approximated  by  27r(0.01) 
by  letting  the  mean  circumference  of  the  annulus  assume  the 
outer  circumference  of  the  unit-disk,  2?r.  The  approximation 
is  99.5%  accurate.  If  the  spatial  correlation  coefficient  has 
a  magnitude  of  90%,  then  the  thickness  of  the  0  annulus  is 
0.1,  the  area  of  the  annulus  can  be  approximated  by  2:r(0.1). 
'Flic  approximation  accuracy  is  95%.  We  see  that  a  combi¬ 
nation  of  the  high  temporal  and  spatial  correlations  yield 
an  rj  of  approximately  (-20  -1-  3)  -f  (- 10  -F  3)  =  -24  dB. 

Next  we  use  Figure  2  to  characterize  some  special  po¬ 
sitions  in  the  isosceles  region.  Hudson’s  textbook[14,  pp. 
52-55]  presented  4  out  of  the  17  characterizations  shown 
there.  These  4  cases  are  denoted  by  the  asterisks. 

The  two  eigenvalues  are  equal  if  and  only  if  what  is  in¬ 
side  the  radical  of  the  eigenvalue  equation  vanishes.  It  can 
be  shown  that  for  this  to  be  true  we  must  have  Tr\/ir2  =  1 
and  that  ^  and  rj  must  be  on  the  left  parabola  of  the 
isosceles-like  region.  The  common  eigenvalue  value  they 
share  is 

A.(^,(l+0')  =  (l+0  =  ^^• 

3. 

Characterizing  Eigenvalues  and  their  Ratio 

We  characterize  the  behavior  of  the  eigenvalues  and 
their  ratio  ^  as  functions  of  the  real  triplet  (^,rj,  ^).  It  is 
useful  to  discuss  the  expressions  for  the  eigenvalue  and  their 
ratio  ^  for  the  special  cases  at  the  apex,  the  baseline,  and 
t!ie  vertical  axis.  But  instead,  we  characterize  the  special 
case  of  equal-strength  arrivals  and  then  the  general  case  will 
be  touched  upon  only  lightly  at  the  end. 

The  small  eigenvalue  is  shown  to  diminish  qualita¬ 
tively  and  quantitatively  for  two-arrival  scenarios  increas- 
in.gly  stressed  with  high  temporal  and/or  spatial  correla¬ 
tions.  The  special  case  of  equal  strength  ^  =  0  dB  arrivals, 
also  important  in  low-angle  radar  tracking,  shares  many  rich 
structures  of  general  It  has  several  additional  unique 
features  for  signal  eigenvalues’  ratio  ^  important  in  array 
processing. 


special  case  of  e<iiial  st,roii(;l.li  =0  111  arrivals 
shares  I  lie  followiiij'  cominoii  feal.iires  of  general  ^  power 
ratios:  straight  lino  contours  for  constant  cigenvahio  A’s  in 
Figure  3,  hyperbolic  X  slices  for  constant  ph:ise-iiulcpenilont 
variable  :y’s  in  f  igure  '1,  parabolic  X  slices  for  constant  phase- 
dependent  variable  ^’s  in  Figure  5. 

However,  we  caution  that  unless  special  elfort  is  made 
in  plotting  results,  the  two  displayed  surfaces  or  curves  may 
not  always  meet  at  the  supposed  places.  This  is  because  the 
boundary  curves  ilefining  parabolas  in  the  f  and  ?/  plane,  in 
general,  do  not  juiss  through  the  grid  points  usetl  in  llic  ( 
and  Tj  plane. 

The  plots  of  the  same  information  or  function,  using 
direct  scale  and  dfl  scale  in  the  vertical  axis  have  their  indi¬ 
vidual  merits.  The  direct-scale  versions  show  the  parabolic 
and  hyperbolic  sections  as  well  as  the  straight  line  contours 
in  their  natural  coordinates  as  derived  from  the  analytical 
studies.  Versions  using  the  dll  scale  better  illustrate  the 
multiplicative  dependence  of  the  eigenvalues  on  some  of  the 
three  independent  parameters.  In  addition,  the  effect  of  the 
smallest  dominant  eigenvalue  is  relative  to  the  threshold  set¬ 
ting,  which  is  frequently  expressed  in  dll  scale;  the  condition 
number  can  be  exhibited  more  compactly  in  this  way. 

As  t;  approaches  0,  the  small  eigenvalue  X2  and  the 
eigenvalue  ratio  ^  tend  toward  -  and  -f  00  dB,  respec¬ 
tively.  Similarly,  as  tj  approaches  0  for  the  equi-power  case, 
^  =  0  dn,  both  the  upper  Ay  and  lower  A2  hyperbolas  tend 
toward  their  asymptotes  intersecting  at  —00  dll.  Since  we 
cannot  display  ±00  dll,  we  choose  to  stop  at  a  small  value 
of  ry  =  — 13  dn,  which  corresponds  to  spatially  orthogonal 
arrivals  with  temporal  correlation  of  97.5%,  or  temporally 
uncorrelated  arrivals  with  a  ^  fractional  beamwidth  spacing 
for  a  (iniform  linear  array. 

We  first  display  for  power  ratio  ^  =  0  dll  the  slices 
of  the  two  signal  eigenvalues  Ai  and  A2  in  dB  for  constant 
phase-independent  variable  rj  in  Figure  6,  and  for  constant 
pliaso-dcpendent  varialilc  ^  in  Figure  7.  While  the  parabola 
sections  are  all  translation  cojiies  of  each  other  for  given 
values  in  the  direct  scale  coordinates,  their  images  in  the  dB 
scale  are  qualitatively  different  because  of  the  nonlinear  na¬ 
ture  of  the  logarithmic  function,  especially  near  zero,  which 
is  so  important  for  our  study.  Next  their  ratio,  which  is  the 
difference  in  dB  of  the  large  and  small  signal  eigenvalues, 
are  displayed  in  Figures  8  and  9. 

We  expect  the  slices  of  the  eigenvalue  ratio  ^  surface 
for  all  constant  tj,  phase-independent  variable,  ancf  for  neg¬ 
ative  constant  phase-dependent  variable,  to  intersect  the 
0-dB  abscissa  axis.  Again,  they  do  not  appear  to  be  so  only 
because  of  the  coarse  grid  systems  used  in  making  these 
plots,  i.e.,  the  left  boundary  curve  of  the  isosceles  region 
does  not  fall  on  the  grid  points. 

Mesh  plots  and  parabolic  contours  for  eigenvalue  ra¬ 
tio  are  shown  in  Figures  10  and  11  respectively.  The 
contours  of  the  eigenvalue  ratio  ^  =  /  form  a  one- 
parameter  family  of  parabolas  having  their  vertices  colo¬ 
cated  at  i(,Tj)  =  (—1,0),  i.e.,  at  the  lower  left  corner  of  the 
isosceles-like  triangle  and  their  common  tangent  the  base¬ 
line  ty  =  0.  This  family  includes  the  left  boundary  of  the 
isosceles-like  triangle  and  the  baseline.  Because  the  contours 
of  the  eigenvalue  ratio  ^  P^^ss  through  the  lower  left 


corner  of  the  iso.sr('leH  like  triangle,  this  comer  is  a  point  of 
disronl  iiiuity  for  I  he  eigenv.ilue  ra(  io  ^  for  t  he  special  case 
oferiual  strength  ^  0  dll.  For  example,  the  two  eigenval- 

ucs  arc  o<]nal  over  ilic  Icfl  parabolic  houiuiary  of  tlic  isoscc- 
les  right-1  riangle-like  region.  I'urt.hcrmore,  the  two  equally 
dominant  eigenvalues  go  down  to  zero  jointly  so  that  the 
condition  number  stays  ns  unity  throughout.  That  is,  the 
condition  numbers  are  the  lowest  possible  there.  However, 
both  eigenvalues  tend  to  zero  as  tj  tends  to  zero,  i.e.,  the 
baseline  of  the  region.  The  limit  of  ^  approaching  from  the 
baseline  of  the  isosceles  region  is  00,  while  that  approaching 
from  the  left  boundary  curve  of  the  isosceles  region  is  0. 

For  the  equi-powered  case,  the  destructive  interference 
scenario  point  is  at  the  left  apex  of  the  isosceles  region.  Both 
eigenvalues  are  zero  and  the  eigenvalue  ratio  ^  nt  this  poii\t 
is  of  the  0/0  form.  'I’his  point  corresponds  to  the  total  can¬ 
cellation  of  signals  at  the  sensor  clement  level.  The  total 
cancellation  happens  not  only  at  the  center  of  the  array  but 
at  every  sensor.  The  two  steering  vectors  co'"'-id'’  and  the 
temporal  waveforms  negate  each  other  completely.  Con¬ 
sider  two  strong  sources  behaving  this  way.  As  the  power  is 
measured  by  turning  off  one  source  at  a  time,  so  the  eq\ial 
power  are  considered  as  high  and  C1U,B  and  asymptotic 
results  apply  when  noise  is  taken  into  consideration. 

We  focus  attention  to  the  condition  number’s  behavior 
over  the  vertical  axis  ^  =  0  of  the  isosceles  region.  For  the 
special  case  of  equal  strength  ^  =  0  dB,  there  is  an  extra 
fi  (IB  for  signal  eigenvalues’  ratio  ^  over  the  vertical  axis 
$  =  0  accounting  for  the  iiresence  o^  an  infinite  slope  of  A’s 
with  respect  to  ry  at  the  apex  for  the  equal-power  arrival 
case.  I'his  is  reflected  in  the  following  two  approximations. 
When  >  10  dB,  we  have 

^(in  (IB)  ss  (  — )(in  dB)  -  ry(in  dB). 

A2 

When  ^  ss  0  dB  and  ry  is  less  than  0.10  (or  -10  dB),  we 
have 

--^(in  (IB)  Ks  — ry(in  dB)  -t-  6. 

A2 

See  Figures  12  and  13. 

The  variation  of  signal  eigenvalue  ratio  ^  across  the 
phase-dependent  variable  ^  is  not  significant.  For  the  special 
ca.se  of  equal  strength  ^  =  0  dB,  there  is  a  d-time  or  6  dB 
increa.se  on  the  eigenvalue  ratio  ^  for  highly  correlated  equi- 
power  arrivals  both  temporally  and  spatially  from  changing 
the  pha-sc  of  the  source  correlation  or  the  angle  difference 
between  the  two  unit-disk  vectors  from  90®  to  0®. 

As  power  ratio  ^  increases  from  unit  value,  the  vari¬ 
ation  of  the  large  eigenvalue  A]  is  essentially  along  the  di¬ 
rection  of  the  horizontal  phase-dependent  variable  The 
smaller  eigenvalue  A2  is  independent  of  both  the  horizon¬ 
tal  phase-dependent  variable  ^  and  ^  near  the  baseline 
for  large  ^  but  is  a  strong  function  o^  the  vertical  phase- 
independent  variable  tj.  The  effect  of  ^  is  essentially  only 
felt  by  the  large  eigenvalue  which  is  not  sensitive  to  whether 
the  scenario’s  (f,7y)  coordinate  is  close  to  the  baseline,  i.e., 
whether  the  arrivals  are  correlated  or  close  to  each  other  as 
far  as  its  order  of  magnitude  is  concerned. 

We  have  used  these  residts  for  assessment  of  scenarios 
in  [16]  for  which  this  paper  is  a  synopsis. 


4.  Conclusion 

Tlic  main  contribution  of  tliis  work  is  a  manageable 
presentation  of  a  compact  map  sliowing  tlic  three  functions, 
Ai,  Aj,  and  over  all  possible  scenarios.  This  enables  one 
to  see  the  relative  positions  among  dillcrenl  scenarios.  We 
also  present  some  easy-to-remember  formulas  that  enable 
one  to  exercise  back-of-envelope  assessment  of  scenarios. 
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I'ignrc  1.  Tlic  nmny-lo-onc- 
ni.'ippliig  Iroin  llic  complex 
p.air  (p,  i)i)  to  real  pair  (£,  t|). 
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I'lgiiio  2.  Mapping  Imm  eaile.sian  piorlncl  of 
Iwo-nnil  rllsk.s  lo  a  pl;mai  legion. 


